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^ C"| Abstract 

Oh In this paper we discuss the law of variation of scale factor a = (t k e l ) « which yields a time-dependent 

deceleration parameter (DP) representing a new class of models that generate a transition of universe from 
the early decelerated phase to the recent accelerating phase. Exact solutions of Einstein's modified field 
equations with perfect fluid and heat conduction are obtained within the framework of Sdez-Ballester scalar- 
tensor theory of gravitation and the model is found to be in good agreement with recent observations. We 
find, for n = 3, k = 1, the present value of DP in derived model as qo = —0.67 which is very near to the 
observed value of DP at present epoch. We find that the time-dependent DP is sensible for the present day 
Universe and give an earmark description of evolution of universe. Some physical and geometric properties 
. of the models are also discussed. 
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1 Introduction 

ly-j | In the last few decades several new theories of gravitation, carefully weighed to be alternative to Einstein's 
theory of gravitation, have been developed according to an orderly plan. In alternative theories of gravitation, 
scalar tensor theories proposed by Brans and Dicke PQ, Nordvedt [2], Wagoner [3], Rose [J, Dun [5], Sdez and 
^vq • Ballester [BJ, Barber [7J, Lau and Prokhovnik [5] are most important among them. There are two categories 
t— I ' of gravitational theories involving a classical scalar field 4>. In first category the scalar field <f> has the dimen- 
sion of the inverse of the gravitational constant G among which the Brans-Decke theory [I] is of considerable 
• i-H ' importance and the role of the scalar field is confined to its effect on gravitational field equations. Brans and 
Decke formulated a scalar-tensor theory of gravitation which introduces an additional scalar field <f> besides the 
metric tensor and a dimensionless coupling constant uj. This theory goes to general relativity for large values 
of the coupling constant uj > 500. In the second category of theories involve a dimensionless scalar field. Sdez 
and Ballester [6J developed a scalar-tensor theory in which the metric is coupled with a dimensionless scalar 
field in a simple manner. This coupling gives a satisfactory description of the weak fields. In spite of the di- 
mensionless character of the scalar field, an anti-gravity regime appears. This theory suggests a possible way to 
solve the missing-matter problem in non-flat FRW cosmologies. The Scalar- Tensor theories of gravitation play 
an important role to remove the graceful exit problem in the inflation era [5]. In earlier literature, cosmological 
models within the framework of Sdez-Ballester scalar-tensor theory of gravitation, have been studied by Singh 
and Agrawal [TUl QT] , Ram and Tiwari [T^] , Singh and Ram [T3] . Mohanty and Sahu [HI [TS] have studied 
Bianchi type-VIo and Bianchi type-I models in Saez-Ballester theory. In recent years, Tripathi et al. [TBJ, Reddy 
et al. [TO EE], Reddy and Naidu [H], Rao et al. [20-22], Adhav et al. |23J, Katore et al. Sahu [25], 

Singh [5BJ , Pradhan and Singh [27] , Socorro and Sabido [25] and Jamil et al. [55] have obtained the solutions in 
Sdez-Ballester scalar-tensor theory of gravitation in different context. Recently, Naidu et al. [30] [31] and Reddy 
et al. [32] have studied LRS Bianchi type-II models in Sdez and Ballester scalar tensor theory of gravitation in 
different context. 
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Recently, Ram et al. [33] obtained Bianchi type-V cosmological models with perfect fluid and heat flow in 
Sdez and Ballester theory by considering a variation law for Hubble's parameter with average scale factor which 
yields constant value of the deceleration parameter. In literature it is common to use a constant deceleration 
parameter [34-40] as it duly gives a power law for metric function or corresponding quantity. But it is worth 
mentioned here that the universe is accelerated expansion at present as observed in recent observations of Type 
la supernova [41-45] and CMB anisotropics [46-48] and decelerated expansion in the past. Also, the transition 
redshift from deceleration expansion to accelerated expansion is about 0.5. Now for a Universe which was de- 
celerating in past and accelerating at the present time, the DP must show signature flipping [49-51]. So, in 
general, the DP is not a constant but time variable. Recently, Pradhan et al. [52] [53] investigated some new 
exact Bianchi type-I cosmological models in scalar-tensor theory of gravitation with time dependent deceleration 
parameter. 

Motivated by the above discussions and observational facts, in this paper, we propose to study Bianchi type- 
V universe with perfect fluid and heat flow in Sdez-Ballcster scalar-tensor theory of gravitation by considering a 
law of variation of scale factor as increasing function of time which yields a time dependent DP. The out line of 
the paper is as follows: In Sect. 2, the metric and basic equations are described. Section 3 deals with the field 
equations and their quadrature solutions. The law of variation of scale factor is given in Sect. 4. Subsection 4.1 
deals with the physical and geometric properties of the universe. Finally, conclusions are summarized in the last 
Sect. 5. 

2 The Metric and Basic Equations 

We consider anisotropic Bianchi type-V line element, given by 

ds 2 = dt 2 - A 2 (t)dx 2 - e 2mx [B 2 (t)dy 2 + C 2 (t)dz 2 ] , (1) 
where A, B and C are metric functions and m is a constant. 

We define the following parameters to be used in solving Einstein's field equations for the metric (JT]). 
The average scale factor a of Bianchi type-V model (Q]) is defined as 

a = {ABC)i. (2) 

A volume scale factor V is given by 

V = a 3 = ABC. (3) 
In analogy with FRW universe, we also define the generalized Hubble parameter H as 

H=- = l(H 1 +H 2 + H 3 ), (4) 
a 6 

where Hi = j|, H 2 = ■§ and H 3 = ^ are directional Hubble factors in the directions of x-, y- and z-axes 
respectively. Here, and also in what follows, a dot indicates ordinary differentiation with respect to t. 

Further, the deceleration parameter q is given by 

ad , , 

We introduce the kinematical quantities such as expansion scalar (9), shear scalar (a 2 ) and anisotropy parameter 
(A m ), defined as follows: 

= (6) 
o 2 = \a^\ (7) 

A m = \±(^)\ (8) 
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where u % = (0, 0, 0, 1) is the matter 4- velocity vector and 
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<Tij = 2 {Ui, a Pf + u jia P?) - -9Pi 
Here the projection tensor has the form 

Pij Qij UiUj. 

These dynamical scalars, in Bianchi type-V, have the forms 

ABC 



9 = 3H 



A B C' 



2a 2 = 



' ■ \ 2 / . \ 2 / ■ \ 2 

A\ ( B\ I C\ 



A I \ B 



C 



(9) 
(10) 

(11) 
(12) 



3 Field Equations and their Quadrature Solutions 

The scalar tensor theories are the generalization of Einstein's theory of gravitation in which the metric is 
generated by a scalar gravitational field together with non-gravitational field (matter). We, consider the simple 
case of a homogeneous but anisotropic Bianchi type-I model with matter term with a scalar field <f>. Our model 
is based on a non-standard scalar-tensor theory, defined in Sdez and Ballester [6] with a dimensionless scalar 
field <p and tensor field g^ . This alternative theory of gravitation is combined scalar and tensor fields in which 
the metric is coupled with a dimensionless scalar field. We assume the Lagrangian 

L = R-u4> k ^' i , (13) 

R being the scalar curvature, cj> a dimensionless scalar field, uj and k arbitrary dimensionless constants and (j> ,% 
the contraction 4>,a9 a ' '■ Here a comma (,) and a semicolon (;) stand for partial and covariant derivative with 
respect to cosmic time t respectively. 



From the above Lagrangian we can establish the action 

I=( (L + 8irL m )(-g)^dx 1 dx 2 dx 3 dx' i , (14) 

where L m is the matter Lagrangian, g is the determinant of the matrix gty, x l are the coordinates, ^2 is an 
arbitrary region of integration. When k = 0, our model is just the Einstein gravity with a massless minimally 
coupled scalar field coupled to gravity. By considering arbitrary independent variations of the metric and the 
scalar field vanishing at the boundary of J^, the variation principle 

51 = 0, (15) 

leads to a generalized Einstein equation 

- u)(j) k (<f>,i(f),j - 7^9ij4>,i<t>' l ' S \ = -87tT-j, 

20V;{ + k^- 1 ^ 1 = 0, (16) 
where Gij = Rij — ^Rgtj is the Einstein tensor; Ty is the stress-energy tensor of the matter Lagrangian L m . 



Since the action / is a scalar, it can be easily proved that the equation of motion 

T% = 0, (17) 

are consequences of the field equations. 
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The energy-momentum tensor is the source of gravitational field through which the effect of the perfect fluid 
with heat flow in the evolution of the universe is performed. The energy-momentum tensor of a perfect fluid 
with heat flow has the form 

Tij = (p + p)uiUj - pgij + hiUj + hjiii, (18) 

where p is the energy density, p the thermodynamic pressure, Ui the four-velocity of the fluid and hi is the heat 
flow vector satisfying 

gnuW = 1, (19) 

and 

h l Ui = 0. (20) 

We assume that the heat flow is in x - direction only so that hi = (hi, 0, 0, 0), hi being a function of time. For 
the energy- momentum tensor (|18[) and Bianchi type-V space-time (fTJ), the Einstein's modified field equations 
(|T6|) . yield the following six independent equations as 

= - p +lcjct> r j> 2 , (21) 

= -p+\u4> r tf 1 (22) 
= -p+JW^, (23) 
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(24) 



AB AC BC A 2 r 2""' 

). (26) 



The law of energy-conservation equation T lJ - — gives 



, A B C\ 2m, 
(P + P)(A + B + c) = A* h - (27) 



Equations (l2"T1) - (f2"4"|) can be written in terms of H, q, a 2 and <f> as 

p = H\2q-l)-a 2 + 1 ^+ l -^^\ (28) 

^3ff 2 -<7 2 -^ + ^ 2 . (29) 

To solve the field equations (l2"Tj) - (l24l) . we follow the well established method of quadrature. For this, subtract- 
ing Eq. from Eq. ^TJ) from (J53J) and Eq. from ([2"3"]). we get the following relations respectively: 

diexp [ki [ ^ ] , (30) 



B 



A . / f dt 



c = d 2 exp yk 2 J ^3 J , (31) 

^ = d 3 exp(k 3 J (32) 

where di,d 2 ,d 3 and ki,k 2 ,k 3 are constants of integration. From Eqs.(|30 p - (|32[) . the metric functions can be 
obtained explicitly as 

A{t)=havcs>(jY /J). (33) 
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/ x , fX 2 f dt\ 
B(t) = l 2 aexp I — s), 



(34) 
(35) 



where 



h = v 7 ^, h = y/d^ds, h = V(d 2 d 3 )-\ 
X 1 = k 1 + k 2 , X 2 = k 3 -h, X 3 = -(k 2 + k 3 ), 
and the constants X%, X 2 , X 3 and lx, l 2 ,l 3 satisfy the relations 
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Figure 1: The plot of deceleration parameter q versus t. 



X 1 + X 2 +X 3 = 0, hhh = l. 
The quadrature expression for the dimensionless scalar field function 0, from Eq. (1261) . is found as 



Mr + 2) /"eft 

,3 



2/(r+2) 



(36) 



(37) 



where <f>o is a constant. 



It is clear from Eqs. (f3"3"|) - (j3"T| that once we get the value of the average scale factor a, we can easily calculate 
the metric functions A, B, C and the scalar function <j). In the next section, we are going to assume the value 
of average scale factor a. 



4 Law of Variation of Scale Factor a = {t k e 1 )^ 

Following Yadav [S3] and Pradhan and Amirhashchi [55], we assume that 

a={t k e t )^ 1 



(38) 
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Figure 2: The plot of isotropic pressure p versus t. Here u> — a — /?2 = X\ = m 



where k and n are positive constants. 

From ([5]) and (|38]) . we obtain the time varying deceleration parameter as 

nk 



q = 



{t + kf 



l. 



(39) 



From Eq. (|3"9")l , we observe that q > for t < Vnk — k and q < for t > \fnk — k. It is observed that 
for n > 3 & k — 1, our model is evolving from decelerating phase to accelerating phase. Also, recent ob- 
servations of SNe la, expose that the present universe is accelerating and the value of DP lies to some place 
in the range — 1 < q < 0. It follows that in our derived model, one can choose the value of DP consistent 
with the observation. Figure 1 depicts the deceleration parameter (q) versus time which gives the behaviour of 
q from decelerating to accelerating phase for different values of (n, k) which is consistent with recent observa- 
tions of Type la supernovae (Perlmutter et al. [41j: Riess et al. [42 } I45 j : Tonry et al. |43j: Clocchiatti et al. [44]). 



Using (15"5)) in Eqs. ([3"31 - (|3"5|) . we obtain the following expressions for scale factors: 



A(t) =l 1 (t k e t ) 1/n exp 

B(t) = l 2 (t k e t ) 1/n exp 

C(t) =/ 3 (t fc e*) 1/ "exp 
Hence the geometry of the universe ([T]) is reduced to 
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ds 2 = di 2_ z 2 (t fe et) 2/„ exp ^2Xi 



dt 

(* fe e*) 3 /" 
dt 

(i fe e*) 3 /" 
dt 
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Figure 3: The plot of energy density p versus t. Here u = a = ft = X\ = m = <po = 1. 



mx l!(t*e*) 2 / n exp 
lf(tV) 2 / n exp 



2X 2 



dt 



3 ./ (i fe e*) 3 / 

2,V ; f dt 

fc e*) 3 A 



dy 2 



(43) 



4.1 Some Physical and Geometric Properties of Model 

In derived model (|43| , the present value of DP is estimated as 

qo = -1 



(44) 



where Hq is the present value of Hubble's parameter and to is the age of the universe at present epoch. If we 
set n = 3 and k = 1 in Eq. (|44|) . we obtain go = —0.67. This value is very near to the observed value of DP 
(i.e., go ~ —0.77) at present epoch (see Cunha et al. [56 ). Hence, we restraint n = 3 and k = 1 in the left over 
discussions of the model and graphical display of physical parameters. 



The solution for scalar function <fi, from (|37jl . is obtained as 

>o(r + 2) f dt 



(£fe e t)3/n 



2/(r+2) 



(45) 



By using the values of the metric functions from eqs. (14T)|) - (14"!?1) into Eq. (|25l) . the expression for the heat flow 
function hi is given by 

hi = 7T77Z ,^/„ , (46) 



3(t* e t)3/n' 



where ft = 2X X - JC 2 - JC 3 . 
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Figure 4: The plot of anisotropic parameter A m versus t. Here a = ftz = 1 
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From Eqs. ([28)) and (f29|) the isotropic pressure (p) and the energy density (p), for model (j43|) . are obtained as 

"♦[(■ 

-2Xi 



2k 



3 / fe 



LU0 n - 
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■ exp 



(i 18 J (t k e l f/ r ' 
dt 



3 / 
P= — 1 + 7 



3m 2 
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(t* e t)3/i 



(47) 



I2( t k e ty/n 



exp 



18/ (t k e t ) 6 / n 

2Xi /■ 



( t fc e i)3/n 



(48) 



In view of (|36|) . it is observed that the above set of solutions satisfy the energy conservation equation ([27]) 
identically and hence represent exact solutions of the Einstein's modified field equations (|2~T j) - (f2"6|) . From Eqs. 
(|47|) and I148|) , we observe that isotropic pressure p and the energy density p are always positive and decreasing 
function of time and both approach to zero as t — > oo. Figures 2 and 3 depict p and p, respectively, versus time 
t showing the positive decreasing function of t and approaching to zero at t — > oo. 

The expressions for physical parameters such as spatial volume (V), directional Hubble parameters (Hi, i — 
1, 2, 3), Hubble parameter (H), scalar of expansion (8), shear scalar (a) and the anisotropy parameter (A m ) for 
model (|43p are, respectively, given by 

V = (t k e t )i, (49) 
X ' (50) 



-If _ 

' ~ n V + iy ' 3(i fc e*)3/™' 



(51) 
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Figure 5: The plot of heat flow function hi versus t. Here f}\ 
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(52) 
(53) 



where /3 2 = X\ 



From Eqs. (|4"9")l and (IB"!]) we observe that the spatial volume is zero at t = and the expansion scalar is 
infinite, which show that the universe starts evolving with zero volume at t = which is big bang scenario. 
From Eqs. (|40|) - (l42l) we observe that the spatial scale factors are zero at the initial epoch t = and hence the 
model has a point type singularity [57]). We observe that proper volume increases with time. 

The dynamics of the mean anisotropic parameter depends on the constant 02 = Xf + X\ + X% . From Eq. 
(|53p . we observe that at late time when t — > oo, A m — > 0. Thus, our model has transition from initial anisotropy 
to isotropy at present epoch which is in good harmony with current observations. Figure 4 depicts the variation 
of anisotropic parameter (A m ) versus cosmic time t. From the figure, we observe that A m decreases with time 
and tends to zero as t — > oo. Thus, the observed isotropy of the universe can be achieved in our model at present 
epoch. 

It is important to note here that lim t _j.o y§z) spread out to be constant. Therefore the model of the universe 
goes up homogeneity and matter is dynamically negligible near the origin. This is in good agreement with the 
result already given by Collins |58j . 

The flow of heat along the x-direction was maximum in early universe, and it diminishes as t — > oo. Figure 
5 describe the variation of heat flow versus cosmic time t which shows the nature of h%. From Eqs. (j46[) and 
(|52|) . we also observe that ?^ = constant which shows that shear scalar is proportional to heat conduction. 
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5 Concluding Remarks 



In this paper we have studied a spatially homogeneous and anisotropic Bianchi type-V space-time within the 
framework of the scalar-tensor theory of gravitation proposed by Sdez and Ballester [B] . The field equations have 
been solved exactly with suitable physical assumptions. The solutions satisfy the energy conservation Eq. (|27|i 
identically. Therefore, new, exact and physically viable Bianchi type-V model has been obtained. To find the 
deterministic solution, we have considered scale factor which yields time dependent deceleration parameters. As 
we have already discussed in Introduction that for a Universe which was deceleration in past and accelerating 
at present time, the DP must show signature flipping [49-51] and so there is no scope for a constant DP. The 
main features of the model are as follows: 

• The model is based on exact and new solutions of Einstein's modified field equations for the anisotropic 
Bianchi type-V space time filled with perfect fluid and heat flow. 

• Our special choice of scale factor yields a time dependent deceleration parameter which represents a model 
of the Universe which evolves from decelerating phase to an accelerating phase. This scenario is consistent with 
recent observations (Perlmutter et al. [41]; Riess et al. [42l|45]; Tonry et al. [43] ; Clocchiatti et al. [44]). 

• Our whole discussions have been concentrated by restraining n — 3, k = 1. By this choice, we find the 
present value of deceleration parameter in derived model as qo — —0.67. This value is very near to the observed 
value of DP (i.e., qo « —0.77) at present epoch (see Cunha et al. [56]). 

• For different choice of n and k, we can generate a class of viable cosmological models of the universe in 
Bianchi type-V space-time. For example, if we set n — 2 in Eq. (|38l) . we find a = V t k e l which is used by 
Pradhan and Amirhashchi [55 in studying the accelerating dark energy models in Bianchi type-V space-time 
and Pradhan et al. [53] in studying Bianchi type-I in scalar-tensor theory of gravitation. If we set k = 1, n = 2 
in Eq. (|38|1 . we find a = Vie* which is utilized by Amirhashchi et al. [59] in studying interacting two-fluid 
scenario for dark energy in FRW universe. If we set k = 1, n = 1 in Eq. (|38|). we find a = te which is exercised 
by Pradhan et al. 60 to study the dark energy model in Bianchi type-U/o universe. It is observed that such 
models are also in good harmony with current observations. 

• It has been observed that lim^o turn out to be constant. Thus the model approaches homogeneity 
and matter is dynamically negligible near the origin. 

• We also observe that £y = constant which shows that shear scalar is proportional to heat conduction (i.e. 
a cc hi. 

Thus, the solutions demonstrated in this paper may be useful for better understanding of the evolution of the 
universe in Bianchi type-V space-time within the framework of Saez-Ballester scalar-tensor theory of gravitation. 
The solutions presented here can be one of the potential candidates to describe the observed universe. 
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